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Abstract 

I summarize the status of neutrino oscillations from world neutrino oscillation data 
with date of October 2005. The results of a global analysis within the three- flavour 
framework are presented. Furthermore, a prospect on where we could stand in neutrino 
oscillations in ten years from now is given, based on a simulation of upcoming long- 
baseline accelerator and reactor experiments. 

1 Introduction 

In the last ten years or so we have witnessed huge progress in neutrino osciUation physics. 
The outstanding experimental results lead to quite a clear overall picture of the neutrino 
sector. We know that there are two mass-squared differences separated roughly by a factor 
of 30, and in the lepton mixing matrix there are two large mixing angles, and one mixing 
angle which has to be small. In the first part of this talk I review the present status of 
neutrino oscillations by reporting the results of a global analysis of latest world neutrino 
oscillation data from solar, atmospheric, reactor and accelerator experiments. This analysis 
is performed in the three-flavour framework and represents an update of the work published 
in Refs. [1,2]. 

The recent developments in neutrino oscillations triggered a lot of activity in the com- 
munity, and many new neutrino oscillation experiments are under construction, or under 
active investigation, to address important open questions, such as the value of the small 
mixing angle 613, leptonic CP violation and the type of the neutrino mass hierarchy. In the 
second part of the talk I try to give an outlook, where we could stand in about ten years 
from now. These results are based on a simulation of up-coming long-baseline accelerator 
and reactor experiments, which are expected to deliver physics results within the anticipated 
time scale [3,4]. 

*Talk given at the XXIX International Conference of Theoretical Physics, "Matter To The Deepest: 
Recent Developments In Physics of Fundamental Interactions", 8-14 September 2005, Ustron, Poland. 
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Table 1; Best fit values (bf), la errors, relative accuracies at Icr, and ia allowed ranges of three-flavour 
neutrino oscillation parameters from a combined analysis of global data, updated from Ref. [2] . 

Three-flavour neutrino oscillations are described in general by the two independent mass- 
squared differences Am2i, Am|]^, three mixing angles 6^12, 6*23, 6^13, and one complex phase 5cv- 
Throughout this work 1 will use the standard parameterization for the PMNS lepton mixing 
matrix 

/ 1 \ / C3 e-''^^sri \ ( C12 S12 \ 

= C23 S23 1 -S12 C12 (1) 

V -S23 C23 / V -e''^^si3 ci3 / V 1/ 

with the abbreviations Sjk = sin 9jk, Cjk = cosOjk- The type of the neutrino mass hierarchy 
is determined by the sign of Am|^: Amg^ > corresponds to the normal hierarchy and 
Amg^ < to the inverted one. 

2 Present status of three- flavour neutrino oscillations 

In Tab. ^ 1 summarize the present status of three-flavour neutrino oscillation parameters. 
The numbers are obtained from a global analysis of current oscillation data from solar [5-7], 
atmospheric [8], reactor [9,10], and accelerator [11] data. Details of the analysis can be 
found in Ref. [2] and references therein. In the following 1 give some brief comments on the 
determination of the 'atmospheric' and the 'solar' parameters, and on the bound on ^13. 

The 'atmospheric parameters'. In Fig.^l show the allowed regions for 6'23 and Am|^ from 
separate analyses of Super-K atmospheric neutrino data [8], and data from the K2K long- 
baseline experiment [11]. The latter probes the z/^ disappearance oscillation channel in the 
same region of Am^ as explored by atmospheric neutrinos. The neutrino beam is produced 
at the KEK proton synchrotron, and originally consists of 98% muon neutrinos with a mean 
energy of 1.3 GeV. The z/^ content of the beam is observed at the Super-K detector at a 
distance of 250 km, where 107 events have been detected, whereas 1511}q have been expected 
for no oscillations. Fig. [H illustrates that the neutrino mass-squared difference indicated by 
the z/^ disappearance observed in K2K is in perfect agreement with atmospheric neutrino 
oscillations. Hence, K2K data provide the flrst conflrmation of oscillations with Amg^ from 
a man-made neutrino source. K2K gives a rather weak constraint on the mixing angle due 
to low statistics in the current data sample, and the constraints on sin^ ^23 of Tab. ^ are 
dominated by atmospheric data. Both data sets give a best fit point of 6^23 = 7r/4, i.e. 
maximal mixing. 

In our analysis of atmospheric neutrino data we neglect the small contribution of os- 
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Figure 1: Allowed regions for sin^ 023 and Am|]^ at 90%, 95%, 99%, and Scr C.L. for atmospheric neutrino 
data (contour lines) and the K2K long-baseline experiment (coloured regions). 

cillations with Am^i. Taking into account this sub-leading effect, in Refs. [12, 13] a small 
deviation from maximal mixing was found, due to an excess of sub-GeV e-like events. This 
indication currently is not statistically significant (about 0.5a), and so-far it has not been 
confirmed by a three- flavour analysis of the Super-K collaboration [14]. 

The 'solar parameters '. In Fig. |2l the allowed regions for 612 and Am^i from analyses 
of solar and KamLAND data are shown. Details of our solar neutrino analysis can be 
found in Ref. [1] and references therein. We use the same data as in Ref. [2] from the 
Homestake, SAGE, GNO, and Super-K experiments [5], and the SNO day-night spectra 
from the pure D2O phase [6], but the CC, NC, and ES rates from the SNO salt-phase are 
updated according to the latest 2005 data [7]. For the KamLAND analysis we are using 
the data equally binned in l/Epj- {Epj- is the prompt energy deposited by the positron), and 
we include earth matter effects and flux uncertainties following Ref. [15] (see the appendix 
of Ref. [2] for further details). We observe from the figure a beautiful agreement of solar 
and KamLAND data. Moreover, the complementarity of the two data sets allows a rather 
precise determination of the oscillation parameters: The evidence of spectral distortion in 
KamLAND data provides a strong constraint on Arn^i, and leads to the remarkable precision 
of 4% at la (compare Tab. Q). As visible in the right panel of Fig. |2l alternative solutions 
around Aml-^ ~ 2 x 10"^ eV^ (~ 1.4 x 10~^ eV^), which are still present in the KamLAND- 
only analysis at 99% C.L., are ruled out from the combined KamLAND+solar analysis at 
about 4(7 {5a). In contrast to Amg^, the determination of the mixing angle is dominated 
by solar data. Especially recent results from the SNO experiment provide a strong upper 
bound on sin^ 6*12, excluding maximal mixing at more than 5a. 

The bound on 613 . For the third mixing angle currently only an upper bound exists. This 
bound is dominated by the CHOOZ reactor experiment [10], in combination with the Am|^- 
determination from atmospheric and K2K experiments. However, recent improved data of 
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Figure 2: Left: Allowed regions for sin^ 6112 and Amji at 90%, 95%, 99%, and Scr C.L. for solar neu- 
trino data (contour lines) and the KamLAND reactor experiment (coloured regions). Right: Combined 
solar+KamLAND analysis. Also shown is the allowed region from solar data only (contour lines). 




Figure 3: Left: Ax^ as a function of sin'^ ^13 for solar+KamLAND, CHOOZ+atmospheric+K2K, and 
global data. Right: Allowed regions for sin^ 612 and sin^ 613 at 90%, 95%, 99%, and Scr C.L. for solar neutrino 
data (solid contour lines), KamLAND (dashed contour lines), and solar+KamLAND data (coloured regions). 
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Label L (E^) trun channel 

Conventional beam experiments: 
MINOS [17] 735km 3 GeV Syr 
ICARUS [18] 732km 17 GeV Syr 
OPERA [19] 732km 17 GeV Syr 
Reactor experiments with near and far detectors: 

D-Chooz [20] l.OSkm ~4MeV 3yr z/g ^ z/g 

Reactor-II [21] 1.70km ~4MeV Syr z/g ^ z/g 
Off-axis super-beams: 

T2K [22] 29Skm 0.76 GeV Syr z/^ z/g, z/^ 

NOz/A [23] 812km 2.22 GeV Syr z/^ ^ z/g, z/^ 

Table 2: Summary of upcoming experiments. 



solar and KamLAND experiments lead to a non-negligible contribution of these experiments 
to the global bound, especially for low values of Amg^^ within the present allowed range [1]. 
From the left panel of Fig. El one deduces the following limits at 90% C.L. {3a): 

( 0.029 (0.067) CH00Z+atm+K2K, 
sin^ ^13 < < 0.041 (0.079) solar +KamLAND, (2) 
[ 0.021 (0.046) global data. 

In the right panel of Fig. IHl we illustrate how the combination of solar and KamLAND 
data leads to a non-trivial bound on Oi^. The allowed regions in the plane of sin^ 6'i2 and 
sin^ 9i3 show the complementarity of the two data sets, which follows from the very different 
conversion mechanisms: vacuum oscillations at a baseline of order 180 km for KamLAND, 
and adiabatic MSW conversion inside the sun for solar neutrinos. For a further discussion 
see the appendix of Ref. [2] or Ref. [16]. 

3 Prospects for the coming ten years 

In this section I discuss the potential of long-baseline experiments, from which results are 
to be expected within the coming ten years. The first results will be obtained by the con- 
ventional beam experiments MINOS and the CNGS experiments ICARUS and OPERA. 
Subsequent information might be available from new reactor experiments. We take as ex- 
amples D-Chooz as a first stage experiment, and a generic second-generation experiment 
labeled Reactor-II, which could be realized at sites in Brazil, China, Japan, Taiwan, or USA, 
see Ref. [21] for an overview. Towards the end of the anticipated time scale results from 
super-beam experiments T2K and NOz/A could be available. The main characteristics of 
these experiments are summarized in Tab. |21 For the simulation the GLoBES software [24] 
is used. Technical details and experiment descriptions can be found in Refs. [3,4,2S]. In 
the following I discuss the expected improvement on the leading atmospheric parameters, 
and the sensitivity to di^, the CP phase 6cp and the type of the neutrino mass hierarchy. A 
discussion of prospects to improve the determination of the leading solar parameters can be 
found in Ref. [26]. 
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Figure 4: Prospective relative errors at 2a on |Ato§]^| (left) and sin^ 023 (right) as a function of the true 
value of Ato§]^ and for the true value sin^ 023 = 0.5. The dots with the error bars indicate the present 
accuracy at 2a from atmospheric and K2K data. The gray shaded region is excluded at 3a by present data. 



The 'atmospheric parameters'. In Fig. |3] I show how the upcoming experiments will 
improve the accuracy on lAmg^l and sin^ 923- Already MINOS will decrease significantly the 
error on |Am3^|. The non-trivial constraint on lAmgj^l from the CNGS experiments results 
from the fact that in this analysis also the disappearance channel is included for these 
experiments (see Ref. [4] for details). For not too small values of |Am|]^| T2K will provide a 
determination of order 2% at 2a, which corresponds roughly to an improvement of one order 
of magnitude with respect to the present error. However, from the right panel of Fig. |3] one 
observes that for most values of lAmg^l T2K will improve the accuracy on sin^ 6'23 only by 
a factor of 2 with respect to the present uncertainty. The main reason for this only modest 
improvement comes from the fact that disappearance experiments measure sin^26'23, and a 
small uncertainty on sin^ 2023 translates into a relatively large error for sin^ 623, if O23 is close 
to maximal mixing. 

It is interesting to note that the lower bound on [Amg^l from NOz/A is comparable to the 
one from T2K, however, the upper bound is significantly weaker because of a strong correla- 
tion between [Amg^l and sin^ ^23- Also the sin^ 6'23 measurement of NOz/A is affected by this 
correlation, which gets resolved only in the range |Am|]^| ^ 3 x 10^'^ eV^. This correlation 
appears because the NOz/A detector is optimized for electrons, whereas the atmospheric pa- 
rameters are determined essentially by the z/^ disappearance channel. Let me add that as 
in Ref. [4] we assume here a low-Z-calorimeter. Using the totally active scintillator detector 
(TASD) as proposed in Ref. [23] improves the performance of NOz/A for the atmospheric 
parameters. 

The limit on 6*13. Fig. shows the limits which can be obtained on sin^ 2^13 by the 
experiments under consideration, if no signal for a finite 613 is found. One observes that 
the conventional beams may improve the present limit roughly by factor of 2, D-Chooz by 
a factor of 4, and the super-beams by a factor of 6. An optimized reactor experiment could 
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Figure 5; Sensitivity to sin^ 26*13 at 90% C.L. The left edges of the bars are obtained for the statistics 
Umits only, whereas the right edges are obtained after successively switching on systematics, correlations, and 
degeneracies, i.e., they correspond to the final sensitivity. The gray shaded region is excluded at 90% C.L. by 
present data. For the true values of the oscillation parameters, we use Am§i = +2.0 X 10-3 eV^ sin^ 26*23 = 1, 
Am|i = 7.0 X 10-^ eV^ sin^ 26*12 = 0.8. 



improve the present bound by more than one order of magnitude, reaching a hmit for sin^ 29is 
below 10~^ at 90% C.L. The limits shown in Fig. |S1 depend on the true value of Amg^; in 
general stronger limits can be reached for larger values of Aml^. 

The figure illustrates that the sin^ 29i3 limits from z/^ appearance beam experiments 
are strongly affected by correlations, mainly between sin^2^i3 and 5cp. The bars labeled 
'Degeneracies' originate from the fact that the hierarchy cannot be determined, and the 
sensitivity is quoted for the hierarchy which gives the worse limit (see the appendix of Ref. [4] 
for a detailed discussion of the sensitivity limit). Clearly, the sin^26'i3 limit from reactor 
experiments is completely free from correlations and degeneracies [3,27], since the Ue survival 
probability does not depend on 6cp and 6*23, and the dependence on the solar parameters 
is negligibly small. These experiments are dominated by systematical uncertainties, related 
mainly to the comparison of the near and far detectors. A possibility to reduce significantly 
the impact of these uncertainties has been presented in Ref. [28]. 

Possibilities for large 613. Finally I discuss the potential of the experiments of Tab.|21if ^13 
is not too far from the present upper bound. In Fig. IHlthe allowed regions in the sin^ 26'i3-5cp 
plane are shown assuming a true value of sin^2^i3 = 0.1. For the super-beam experiments 
T2K and NOi/A one observes the strong correlation between sin^ 2^13 and 6cp, which intro- 
duces a large uncertainty in sin^2^i3 but does not permit to draw any conclusions on 6cp- 
In contrast, the reactor experiment provides an accurate measurement of sin^ 2^13 with an 
accuracy at the level of 10% at 90% C.L. None of the experiments on its own can identify 
the mass hierarchy. As indicated in the figure, the best fit point of the wrong hierarchy is as 
good as the true best fit point. However, from the combination of T2K+N0z/A-|- Reactor-II 
the wrong hierarchy can be disfavoured with Ax^ = 3.1, which corresponds roughly to the 
90% C.L. The crucial element for the mass hierarchy determination is the long baseline of 
812 km for NOz/A, leading to matter effects which allow to distinguish between normal and 
inverted hierarchy. Let me add, that the hierarchy sensitivity strongly depends on the true 
value of Sep] typically it is optimal for 6cp — —90°, and worst for 6cp — 90° [4]. 
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Figure 6: The 90% C.L. (solid) and 3a (dashed) allowed regions (2 d.o.f.) in the sin^ 20i3-(5cp plane 
for the true values sin^ 26*13 = 0.1 and ^cp = 90°. The black curves refer to the allowed regions for the 
normal mass hierarchy (assumed to be the true hierarchy), whereas the gray curves refer to the agn{Arn^i)- 
degenerate solution (inverted hierarchy), where the projections of the minima onto the sin^ 2f?i3-(5cp plane 
are shown as diamonds (normal hierarchy) and triangles (inverted hierarchy). For the latter, the Ax^-value 
with respect to the best-fit point is also given. 

As visible from Fig.lHlno information can be obtained on CP violation, i.e., at least one of 
the CP-conserving values 6cp = 0, 180° is contained within the 90% C.L. region, even though 
the assumed true value Sep = 90° corresponds to maximal CP violation. The reason is that 
no anti-neutrino data is included in this analysis, since due to the low cross sections it seems 
unlikely that significant anti-neutrino data will be available within the anticipated time scale. 
Note however, that from the combined analysis some values of 6cp can be excluded for a 
given mass hierarchy. 

4 Conclusions 

I have reviewed the present status of neutrino oscillations from world neutrino oscillation 
data, including solar, atmospheric, reactor and accelerator experiments. The results of a 
global analysis within the three-flavour framework have been presented, and in particular, 
the bound on 613, which emerges from the interplay of various data sets has been discussed. 
Furthermore, a prospect on where we could stand in neutrino oscillations in ten years from 
now has been given. Based on a simulation of upcoming long-baseline accelerator and re- 
actor experiments the improvements on the leading atmospheric parameters, as well as the 
sensitivity to ^13, 6cp and the neutrino mass hierarchy have been discussed. 
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